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The low energy regime of cosmological BPS-brane configurations with a bulk scalar field is stud- 
ied. We construct a systematic method to obtain five-dimensional solutions to the full system of 
equations governing the geometry and dynamics of the bulk. This is done for an arbitrary bulk 
scalar field potential and taking into account the presence of matter on the branes. The method, 
valid in the low energy regime, is a linear expansion of the system about the static vacuum solution. 
Additionally, we develop a four-dimensional effective theory describing the evolution of the system. 
At the lowest order in the expansion, the effective theory is a bi-scalar tensor theory of gravity. 
One of the main features of this theory is that the scalar fields can be stabilized naturally without 
the introduction of additional mechanisms, allowing satisfactory agreement between the model and 
current observational constraints. The special case of the Randall-Sundrum model is discussed. 



I. INTRODUCTION 

Brane-world models have become an important sub- 
ject of research in recent years. The simple and intuitive 
notion that the standard model of physics is embedded in 
a hypersurface (brane), located in a higher dimensional 
space (bulk) , has significant implications for gravity, high 
energy physics and cosmology. It is strongly motivated 
by recent developments in string theory and its extended 
version, M-theory [H, S • Additionallv^it has been pro- 
posed on a phenomenological basis [1, H; 13 ^o address 
many questions which cannot be answered within the 
context of the conventional standard model of physics, 
such as the hierarchyproblem and the cosmological con- 
stant problem 0,|lElISIIllIl- 

The basic idea is that gravity -and possibly other 
forces- propagate freely in the bulk, while the Standard 
Model's fields are confined to a 4-dimensional brane (see 
[13L I14L llSl Tq for comprehensive reviews on brane-world 
models). In this type of model, constraints on the size of 
extra dimensions are much weaker than in Kaluza-Klein 
theories [i3|, although, Newton's law of gravity is still 
sensitive to the presence of extra dimensions. Gravity has 
been tested only to a tenth of millimeter [l^, therefore, 
possible deviations to Newton's law, below that scale, 
can be envisaged. In general, one of the main modifi- 
cations to general relativity offered by brane models, is 
the appearance, in the effective Einstein's equations of 
the system, of novel terms that depend quadratically on 
the energy-momentum tensor of the matter content of 
the observed Universe. When these terms can be ne- 
glected the brane system is in the low energy regime. 
It is well known that the most recent cosmological era 
must be considered within the low energy regime. De- 
spite this fact, there are still important modifications to 
general relativity that must be considered, and current 
observational constraints are highly relevant for the phe- 
nomenology of brane models. 

Cosmologically speaking, the evolution of the branes 



in the bulk is directly related to the evolution of the ob- 
served Universe [l^ |23|. So far, most of cosmological 
considerations of brane-worlds have been centered on the 
models proposed by Randall and Sundrum J3, Il0| • They 
considered warped geometries in which the bulk-space 
is a slice of Anti de Sitter space-time. One of the main 
problems of these scenarios is the stabilization of new ex- 
tra degrees of freedom appearing from the construction. 
More precisely, when the four-dimensional equations of 
motion for the branes are considered, a new scalar de- 
gree of freedom, the radion, emerges. The stabilization 
of the radion is an important phenomenological prob- 
lem, and many different mechanisms to stabilize it have 
been proposed 21, 23|. Clearly, brane models which only 
consider the presence of gravity in the bulk, such as the 
Randall-Sundrum model, are not fully satisfactory from 
a theoretical point of view. There is no reason a pri- 
ori for gravity to be the only force propagating in the 
bulk, and other unmeasured forces may also be present. 
For instance, string theory predicts the existence of the 
dilaton field which acts as a massless partner of the gravi- 
ton. If such forces are present in the bulk, they should be 
suppressed by some stabilization mechanism in order to 
evade observation (similar to the case of the radion field 
in the Randall-Sundrum model). The problem of finding 
a stabilization mechanism is an active field of research. 

In this paper we are going to consider a fairly general 
brane-world model (BPS-branes), which has been mo- 
tivated as a supersymmetric extension of the Randall- 
Sundrum model [M 11 El m 113 (see [M El for 
phenomenological motivations). It consists of a five- 
dimensional bulk space with a scalar field -0, bounded 
by two branes, cri and 02- The main property of this 
system is a special boundary condition that holds be- 
tween the branes and the bulk fields, which allows the 
branes to be located anywhere in the background with- 
out obstruction. In particular, a special relation exists 
between the scalar field potential U{ip), defined in the 
bulk, and the brane tensions C/b(V'^) and Ub{4''^) de- 



fined at the position of the branes. This is the BPS con- 
dition. Due to the complexity of this set up, the main 
efforts of previous works have been focused on the special 
case Up {fp) oc e"'^ , where a is a parameter of the theory 
[sol Isijl . This particular potential is well motivated by 
supergravity in singular spaces and string theory. For 
instance, when a = 0, the Randall-Sundrum model is 
retrieved, meanwhile, for larger values, a = 0(1), one 
obtains the low-energy effective action of heterotic M- 
theory, taking into account the volume of the Calabi-Yau 
manifold 32] . Nevertheless, other classes of potentials 
also have physical motivation and should be considered. 
For instance, in 5-dimensional supergravity, one can ex- 
pect the hnear combination Ub{^) = Vie"^"^ + ¥26"^"^, 
where ai and 02 are given numbers |2,'t| . It must be 
stated, however, that up to now the problem of assuming 
an arbitrary potential Usifp) has not been considered in 
detail, and therefore, the low energy regime of this type 
of system is not well understood. For example, it is not 
known how the gravitational fields nor the bulk scalar 
field behave near the branes, and the general belief is 
that these issues should be investigated via numerical 
techniques. 

Here we study the low energy regime of brane-models 
with a bulk scalar field. One of the most important as- 
pects is that we shall not restrict our treatment to any 
particular choice of the BPS potential, allowing a clearer 
understanding of 5-dimcnsional BPS systems, and focus 
our study on the problems mentioned previously. In par- 
ticular we shall develop a systematic method to obtain 
5-dimensional solutions for the gravitational fields and 
the bulk scalar fields. Additionally, by adopting the pro- 
jective approach, we develop an effective 4-dimensional 
theory, valid in the low energy regime, and indicate that 
it is equivalent to the moduli-space approximation. This 
effective theory is a bi-scalar tensor theory of gravity of 
the form: 




^g'^-iabd^rd.^p' - ^V 



+Si[^i, Afg^,] + S2[<5^2, Alg^,], 



(1) 



where a = 1, 2 labels the branes, ip'^ corresponds to the 
projected bulk scalar field at the brane aa, Jab is a a- 
model metric of the bi-scalar theory, 5*1 and S2 are the 
actions for matter fields ^1 and ^2 living in the respec- 
tive branes, and Ai and A2 are warp factors dependent 
on the scalar fields. Obtaining this effective theory is 
a great achievement that allows the study of this class 
of models within the approach and usual techniques of 
multi-scalar tensor theories [33l l34l| . 

This paper is organized as follows: In Section ^ we 
review brane systems with a bulk scalar field. There, we 
illustrate in a more technical way some of the problems 
that this paper addresses. In Section IIIII we work out 
and analyze the static solutions for the vacuum configu- 
ration of BPS-branes. This enables us to construct the 
low energy regime expansion in Section IIVI Then, in 
Sectional we develop the four-dimensional effective the- 



ory valid in the low energy regime. There, we will focus 
our attention on the zeroth order effective theory which 
is a bi-scalar tensor theory of gravity. In Section IVII 
we study the cosmology of the effective theory. For in- 
stance, we consider the constraints coming from current 
observations on deviations to General Relativity and an- 
alyze the stabilization problem of BPS-systems. It will 
be found that, in order to stabilize this type of system 
it is enough to have the appropriate potential Ub, and 
therefore no exotic mechanisms are required in general. 
The conclusions are summarized in Section |VII| 



II. 



BPS-BRANES WITH A BULK SCALAR 
FIELD 



We start this section with the introduction of 
5-dimensional BPS configurations, deriving the 5- 
dimensional equations of motion governing the dynamics 
of the system and deducing the projected equations at 
the position of the branes. These equations depend on 
the configuration of the bulk as well as on the matter con- 
tent in the branes; we analyze the difficulties arising from 
this dependence to find solutions to realistic cosmological 
configurations. 



A. Basic Configuration 

Let us consider a 5-dimensional manifold M provided 
with a coordinate system X^, with A ~ 0, ... ,4. We 
shall assume the special topology AI — a x S^ /Z2, where 
cr is a fixed 4-dimensional lorentzian manifold without 
boundaries and S^ IZ2 is the orbifold constructed from 
the 1-dimensional circle with points identified through a 
Z2-symmetry. The manifold M is bounded by two branes 
located at the fixed points of S^ IZ2. Let us denote the 
brane-surfaces by a\ and 02 respectively. We shall usu- 
ally refer to the space bounded by the branes as the bulk 
space. We will consider the existence of a bulk scalar 
field 1/' living in M with boundary values, 1/'^ and V'^, at 
the branes. Consider also the presence of a bulk poten- 
tial U{^) and brane tensions Vi('0^) and V2(V'^) (which 
are potentials for the boundary values V'^ and V''^)- Ad- 
ditionally, we will consider the existence of matter fields 
^1 and ^2 confined to the branes. Figure ^ shows a 
schematic representation of the present configuration. 

Given the present topology, it is appropriate to intro- 
duce foliations with a coordinate system x^ describing a 
(as well as the branes a\ and 02) where /i = 0, . . . , 3. Ad- 
ditionally, we can introduce a coordinate z describing the 
S^ /Z2 orbifold and parameterizing the foliations. With 
this decomposition the following form of the line element 
can be used to describe M: 



ds' 



{N^ + NPNp)dz^ + 2N^dx''dz + gf.^dx'^dx". (2) 



Here N and N^ are the lapse and shift functions for the 
extra dimensional coordinate z and, therefore, they can 
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FIG. 1: Schematic representation of the 5-dimensional brane 
configuration. In the bulk there is a scalar field tp with a bulk 
potential U{ip). Additionally, the bulk-space is bounded by 
branes, ai and cr2, with tensions given by Vi{tp^) and V2{'4'^) 
respectively, where tp'^ and tp'^ are the boundary values of ip. 



be defined up to a gauge choice (for example, the gaussian 
normal coordinate system is such that N'^ ~ 0). Addi- 
tionally, (7^1/ is the puUback of the induced metric on the 
4-dimensional foliations, with the (— , -|-, -I-, -|-) signature. 
The branes are located at the fixed points of the S^ /Z2 
orbifold, denoted hy z = zi and 2: = 2:2, where 1 and 2 
are the labels for the first and second branes, ci and CT2- 
Without loss of generality, we shall take zi < Z2- 
The total action of the system is 



■^G + "Sbulk + 'S'bR, 



(3) 



where Sq is the action describing the pure gravita- 
tional part and is given by Sg — Seh + Sqh, with 
5eh the Einstein-Hilbert action and Squ the Gibbons- 
Hawking boundary terms; S'buik is the action for the mat- 
ter fields living in the bulk, which can be decomposed 
into S'buik = Sji, + Sm. Here S^ is the action for the bulk 
scalar field ip (including its boundary terms), while Sai 
is the action for other matter fields (in the present work 
we will not specify this part of the theory) . Finally Sbr 
is the action for the matter fields, ^1 and ^2, confined 
to the branes. 

Let us see in more detail the different contributions ap- 
pearing in equation Q. First, note that with the present 
parameterization for M, we can write an infinitesimal 
proper element of volume as dV = dzd^x^—gN 
it is possible to express Sq in the following way 



Then, 




24 



(4) 
-^N{R-[K^,K^^''-K^]), 



where R is the four-dimensional Ricci scalar con- 
structed from 5^,^, and Kg ~ SvrGs, with G5 
the five-dimensional Newton's constant. Additionally, 
K^u — \g'^^ — y ^Ni, — Vi/A^^] /2N is the extrinsic cur- 
vature of the foliations and K its trace (the prime denotes 
derivatives in terms of z, that is ' = dz, and covariant 
derivatives, V^, are constructed from the induced metric 
Qfj^^ in the standard way). On the other hand the action 
for the bulk scalar field, S^, can be written in the form 



S,b = 



^■4>: 



dz d X\ 



~'^j;p df,tp dzip 




j^) +id^r 



c2 



where Si and Si are boundary terms given by 






2«^5 Ja 



d^x^^VM^), 



^^^5 Jan 



(5) 



(6) 



(7) 



at the respective positions, zi and ^2- In the above ex- 
pressions U[tp) is the bulk scalar field potential, while 
Vi{^^) and V2(V'^) are boundary potentials defined at 
the positions of the branes (also referred as the brane 
tensions). In the present case (BPS-configurations), we 
shall consider the following general form for the poten- 
tials: 



U = U + u, 
Vi = Ub + vi, 

V2 = Ub + V2, 



(8) 

(9) 

(10) 



where U and Ub are the bulk and brane superpotentials, 
and the potentials u, vi and V2 are such that \u\ <C \U\ 
and |wi|, 1^2 1 <C \Ub\- In this way, the system is domi- 
nated by the superpotentials U and Ub- The most impor- 
tant characteristic of this class of system is the relation 
between U and Ub (the BPS-relation), given by: 



U^{d.4,UBf -Ul. 



(11) 



As we shall see later in more detail, the BPS property 
mentioned earlier (the fact that the branes can be lo- 
cated anywhere in the background, without obstruction) 
comes from the relation above. This specific configura- 
tion, when the potentials u,vi,V2 =0 and no fields other 
than the bulk scalar field are present, is referred as the 
BPS configuration. It emerges, for instance, in models 
like N = 2 supergravity with vector multiplets in the bulk 
|25|. More generally, it can be regarded as a supersym- 
metric extension of the Randall-Sundrum model. In par- 
ticular when Ub is the constant potential, the Randall- 
Sundrum model is recovered with a bulk cosmological 
constant A5 = (3/8)C/ = — (3/8)C/^, and the usual fine 
tuning condition is given by equation Hll|) . The presence 
of the potentials u, vi and V2 are generally expected from 



supersymmetry breaking effects. However, their specific 
form is not known and, up to now, they must be arbi- 
trarily postulated 35] . 

Finally, for the matter fields confined to the branes, we 
shall consider the standard action: 



Sbr = Si[^i,g^^{zi)] + S2[^2,gfiiyiz2) 



(12) 



where ^i and ^2 denote the respective matter fields, and 
gfii^iza) is the induced metric at position Zq. In this arti- 
cle we shall assume that the matter fields are minimally 
coupled to the bulk scalar field. That is to say, we will 
assume that Sbr does not depend explicitly on ip. Never- 
theless, it is important to point out that, in general, one 
should expect some nontrivial coupling between the bulk 
scalar field -0 and the matter fields ^1 and '^2, that could 
lead to the variation of constants in the brane- world sce- 
nario iSGl. Despite these comments, we show that, in a 
minimally coupled system, an effective coupling is gener- 
ated between the matter fields and the bulk scalar field. 



B. Equations of Motion 

In this subsection we derive the equations of motion 
governing the dynamics of the fields living in the bulk. 
These equations are obtained by varying the action S'tot 
with respect to the gravitational and bulk scalar fields, 
without taking into account the boundary terms. These 
will be considered in the next subsection, in order to 
deduce a set of boundary conditions for the bulk fields. 

The 5-dimensional Einstein's equations of the system 
can be derived through the variation of the total action 
S'tot with respect to N, Nf^ and g^i,. Respectively, these 
equations are found to be: 

R+[K^,,K'^''-K^]=24S, (13) 



V^ m - KSi^] = -nl J„ (14) 



^ LiU 



-^fiiy 



-g^, [Kp^RP^ - if 2] ^ ^^^^^ 



-K'^^K, 



— Kk 

3 " 



^f^i' .9i-iv\^cj 



35) 



(15) 



The above set of equations deserves a detailed descrip- 
tion. G^n , in (|15|) , is the Einstein tensor constructed from 
the induced metric g^^. The quantities T^i/, J^, and S 
are the 4-dimensional stress energy-momentum tensor, 
vector and scalar respectively, characterising the matter 
content in the bulk space. They are defined as 



S = -- 



1 SShu 



Ik 



T = 



V-9 SN 

1 (^Sbulk 
2 (^Sbulk 



(16) 



where Shuik — 8^1, + S^ is the action for the matter fields 
living in the bulk. The above quantities are decomposi- 
tions of the usual 5-dimensional stress energy-momentum 
tensor. Tab- To be more specific, the definitions in l|l()|) 
satisfy: T^.^X^^X^Tab, J^. = -X^^u^Tab. and 
S = —n^n^ Tab, where X^ is the puUback between 
the coordinate system x^ of the foliations and the 5- 
dimensional coordinate system X^ of M, and where n^ 
is the unit-vector field normal to the foliations. In par- 
ticular, the contributions coming from the scalar field 
action, S^, are given by: 



^(■0) 



Ji'^) = 






1 N^^ 



8«i 



44 



1 / 1 /VM 



+ idi>Y+U] -dpi;d,^ 



(17) 



Also in H15|) is S^t/, the projection of the 5-dimensional 
Weyl tensor, C^adb^ ^^ the foliations given by E^^^ = 
X^ X^^ nc n^ d adb- ^^ ^^'^ present set up, this can be 



n as: 






E^v - 


-- ir^if„, - 1 (V^V.TV + £^ 




3 ^ 


f^,~\g^,{f:+~S) 



^K^,,) 



(18) 



where £Ar„ is the Lie derivative along the vector field 
Nn^. At this point it is convenient to adopt the gauge 
choice iV^ = 0, which corresponds to the choice of 
gaussian normal coordinates. With this gauge, one has 
£Nn = dz, and the treatment of the entire system is 
greatly simplified. Next, we can obtain the equation of 
motion for the bulk scalar field ip. Varying the action 
with respect to ■0 gives: 



g^'-V^ (Nd^^P) + i^'/N)' - 



-*'4S- <-) 



(Recall that we are now using N^ = 0). 

It is possible to see that equations ((T^ and ((T5|l are 
the dynamical equations for the gravitational sector (no- 
tice that second order derivatives in terms of z are only 
present in the projected Weyl tensor E'^^). Meanwhile, 
equation (|14|l constitutes a restriction for the extrinsic 
curvature K^^, in terms of the matter content of the bulk. 
As we shall soon see in more detail, the conservation of 
energy-momentum on the branes is a direct consequence 
of equation ifT^ . 



C. Boundary Conditions 

We can now turn to the boundaries of the system. The 
former variations leading to equations H13|l - H15|) and H19|) 



also give rise to boundary conditions that must be re- 
spected by the bulk fields at positions Zi and Z2- They 
emerge when we take into account the boundary terms 
present in the action (including the matter fields in the 
branes). For the first brane, cti, these conditions are 
found to be: 

{K^, - g^.KJU, - 1[Ub+ t'l] 5m- - ^^T^^, (20) 

lW/N}\zi - d^pB+vi], (21) 

while for the second brane, a2, these conditions are: 






-i^[UB+V2]g^u-4TU, (22) 



-d^ [Ub + V2] 



(23) 



Equations (|20|l and H22|) are the Israel matching condi- 
tions and equations (|21l) and H23|) are the BPS match- 
ing conditions for the bulk scalar field. In the previous 
expressions the brackets denote the difference between 
the evaluation of quantities at both sides of the branes, 
i-e {f}\z = lim£^o[/(2 + e) - I{z - e)]. Additionally, 
T^j, denotes the energy-momentum tensors of the matter 
content in the brane da, and is defined by the standard 
expression: 



rpa 



SSa 



^Sgf^-' 



(24) 



Since we are considering an orbifold with a Z2-symmetry, 
and the branes are positioned at the fixed points, the 
conditions H2U|) - (|23|I can be rewritten as 

i^M^' - 5m-^ - i\Ub+ t'l] 5p. - ^Ti , (25) 



ij' = Nd^pB+vi], 



(26) 



for the first brane, cti, and 



^M- - g^..K - j[Ub+ V2] 5p. + -i^T^,, (27) 



iP' = Nd4,[UB+V2], 



(28) 



for the second brane, (72. These conditions are of utmost 
importance. They relate the geometry of the bulk with 
the matter content in the branes. In the next subsection 
we shall examine in more detail how these conditions are 
related to the 4-dimensional effective equations in the 
branes and analyze some of the open questions. 



D. Short Analysis and Open Questions 

Equations (|13|I - H15(I and 1)1911 can be evaluated at the 
position of the branes, with the help of the boundary 
conditions H25|l - H28|) . to obtain a set of 4-dimensional ef- 
fective equations. This is the projective approach. For 



instance, evaluating equations H13|l and H15() at the posi- 
tion of the first brane, we find: 



R 



'-ViT - g^'^Hp 



-W 



-id^f 



Gfiv — —r^iTfiu 






e 



/j-i> 



E, 



d^il^d^ijj - -g^^^idipY' 



3 

16 



(29) 



g^uW, (30) 



where 0^i/, II^j, and W have been defined as 
2 



0(m) 






\9,. {f^n^ 



35(" 



n 



fJ.V 



1, 



-TT.^^T. 



rp(T 



' ^giiv^cp^ „.9/Ji'-' 



w = u 



^1 -.'■ 



(31) 

(32) 
(33) 



Notice that we have dropped the index labelling the first 
brane in the energy momentum tensor (that is, we have 
taken T^^ = T^^ and T = g'^'T^^). Similarly, the pro- 
jection of equation H19|l on the first brane leads us to the 
4-dimensional effective bulk scalar field equation. This is 
found to be: 



D^j 



1 

'n 



d^'NdJj 



6 9^ 



T- Ai/;- 



ldu_ 






■Vi 



>' 



d^^ 



(34) 



where we have defined the loss parameter Aip of the sys- 
tem as: 



1 



Ai; = -[i;'/N^d.^UB]'. 



N 



(35) 



Additionally, from the evaluation of equation 1)141) at po- 
sition zi, one finds the energy-momentum conservation 
relation for the matter fields of the first brane: 



v„r.^ = 2j('"). 



(36) 



Note that if u = ui = then W = 0. The precise form 
of f^^\ Ji"^ and ^("), in equations (P ^ . P7 |) and P^ . 
depend on the bulk matter fields considered in Sm- From 
now on we will not consider any contribution from Sm 
and therefore we shall take f^'°^ = J^""^ = ^('") = (that 
is, the only fields living in the bulk are the gravitational 
fields (7^1, and N, and the bulk scalar field ip). 

Equations (P5|l . (|^ and (|S^ are the projected equa- 
tions describing the theory in the first brane (an analog 
set of equations can be obtained for the second brane). 
To solve them, it is important to know the precise form 



of the projected Weyl tensor _Ep^ and the loss parameter 
A^. In order to find a complete solution of the entire 
system it is necessary to find a solution for these quan- 
tities. They contain the necessary information from the 
bulk, or equivalently, from a brane frame point of view; 
they propagate the information from one brane to the 
other. A consistent method to compute E^^, and Aijj is 
an open question. 

Observe that the 4-dimensional effective Newton's con- 
stant, G4, can be identified as 87rG4 = K5V1/4 and, there- 
fore, it is a function of the projected bulk scalar field i/"- 
Additionally, it can be seen that there are two scalar de- 
grees of freedom. Naturally, one them is the projected 
bulk scalar field V'- The second one corresponds to the 
projected value of lapse function N, which is normally 
referred as the radion field. In the absence of supersym- 
metry breaking terms u, vi and V2, these scalar degrees 
of freedom are massless and are the moduli fields. The 
stabilization of moduli fields is one of the most impor- 
tant phenomenological issues in brane models to ensure 
agreement with observations. 

The term 11^,^, in equations H3U|I and (|34|l . contains 
quadratic contributions from the energy momentum ten- 
sor Tfj_,y. Therefore, II^i, is relevant at the high energy 
regime and it will turn out to be important for the early 
universe physics. At more recent times, in order to agree 
with observations, this term must be negligible, which is 
only possible if K5 \Ub\ ^ |2^|. This is the low energy 
regime. Neglecting fl^i, makes the previous equations of 
motion take the form: 



R 



^ liV ^— 



-^UBT+\{d^)\ 



—7- Ub T)iiy — E^j^i, 



d^,^d^ip- -g^„{dTp)' 



(37) 



(38) 



D^ 



±d'^Nd,^-^^T-Ai,^ (39) 

A* 6 oyj 



where we have, for simplicity, also dropped the terms con- 
taining u and vi . Recall the additional equation express- 
ing energy-momentum conservation on the first brane, 
which now reads: 



V^T,^ 



0. 



(40) 



This version for the equations of motion are much simpler 
than the full system, however, the difficulties mentioned 
above are still present and many questions remain to be 
answered. For example, the usual approach is to neglect 
the effects of the loss parameter A^; and take it to be 
zero. If this is the case, observe that the projected scalar 
field ip, in equation (|39|) . is driven by the derivative d,pUB 
times the trace of the energy-momentum tensor. One 
could therefore expect the scalar field to be stabilized 
provided the appropriate form of Ub- This possibility, 
however, will depend on the dynamics of the radion field 
N , which, at the same time depends on the dynamics 



of the inter-brane system. The stabilization of the scalar 
field and the radion are usually solved by the introduction 
of new potential terms for ip at the branes and the bulk. 
In the present case that could be the role of u, vi and V2- 
In the rest of the paper we shall focus our efforts to 
solve these problems in the low energy regime. For in- 
stance, we will see that it is possible to introduce a proce- 
dure to compute E^^, and Aip with any desire precision. 
Additionally, we show that it is possible to stabilize the 
scalar degrees of freedom solely with the introduction of 
the appropriate potential Ub (that is, without the need 
of supersymmetry breaking terms u, vi and W2). 



III. STATIC VACUUM CONFIGURATION 

Before studying realistic solutions in detail, namely 
when matter fields are present in the branes, let us ana- 
lyze the static vacuum solutions to this system. This will 
be highly relevant for the development of the low energy 
expansion and the effective theory, to be worked out in 
the next sections. 



A. Static Vacuum Solution 

Let us consider the presence of the bulk scalar field V' 
and gravitational fields A'' and g^^. In this case, we are 
mainly concerned with static vacuum solutions to the 
system of equations (I13|l - (|15|) and H19|) . Therefore, we 
shall assume that tp, N and g^i, are such that: 



du^ = 0, 



df^N - 



0, and Gn 



0, 



(41) 



where G^^ is the Einstein's tensor constructed from g^^. 
To solve the system, it is sensible to consider the following 
form for the induced metric: 



5m"^ '=^^'^{z)9t^y{x), 



(42) 



where g^i^(x) is a metric that depends only on the space- 
time coordinate a;, and which necessarily satisfies the 4- 
dimensional vacuum Einstein's equations G^^ = 0. In 
this way, all the dependence of the induced metric on 
the extra-dimensional coordinate z is contained in the 
warp factor (jj{z). With this form of the metric, the set 
of matching conditions (|25|I - H28(I can be rewritten in the 
following simple way: 






1 



A^ 



^NUb. 
OUb 



9V 



(43) 
(44) 



More importantly, it is possible to show that these rela- 
tions solve the entire system of equations (tT^ - ((T3|) and 
(|19|1 . This important fact constitutes one of the main 
properties of BPS-systems. It means that the branes can 
be arbitrarily located anywhere in the background, with- 
out obstruction. It should be clear that when matter is 



allowed to exist in the branes, the boundary conditions 
(|43|l and H44I) will not continue being solutions to the sys- 
tem of equations l|13|l - H15|l . and the static configurations 
will not be possible in general. The introduction of mat- 
ter in the branes will drive the branes to a cosmological 
evolution. 

In the static vacuum solution expressed through equa- 
tions (|43|) and H44|) , the dependence of the lapse function 
iV, in terms of the extra dimensional coordinate z, is 
completely arbitrary, though it must be restricted to be 
positive in the entire bulk, and its precise form will corre- 
spond to a gauge choice. Let us assume that ^(z) satisfies 
equations (|43ll and (|44ll . and that it has boundary values 
ip^ and ip'^ defined by: 



i''=^{zi) 



and 



V-' = ^(22). 



(45) 



Since we are interested in the static vacuum solution, ip'^ 
and tp'^ are just constants. The precise form of "0(2:), as 
a function of z, depends on the form of C/b('0) and the 
gauge choice for TV. However, it is not difficult to see that 
tp'^ and "0^ are the only degrees of freedom necessary to 
specify the BPS state of the system. That is, given a 
gauge choice for N, in general we have: 

^^ip{z,^\ip'^), and N ^ N{z,ip\ip'^). (46) 

Observe, additionally, that u!{z) can be expressed in 
terms oi ip{z) in a gauge independent way. Using both 
relations, (|43|) and H44|) . we find: 



oj{z) = exp 



1 



'AW 



1/.1 



a ^{tp)d^p 



M) = 



1 OUb 

Ub dtp ' 



(47) 
(48) 



In the last equations we have normalized the solution 
uj{z) in such a way that the induced metric to the first 
brane is g^i^, though this choice is not strictly necessary. 
The induced metric on the second brane is, therefore, 
conformally related to the first brane, with a warp factor 

UJ{Z2)- 

Let us emphasize the fact that, given a solution 'ip{z) 
and a gauge choice N, satisfying equations H43|l and (|44|l . 
the entire system can be specified by providing the de- 
grees of freedom ip^, ip"^ and g^^. In the next section 
we shall see that when matter fields are considered, the 
vacuum solution is perturbed and ip^ , ip'^ and g^i, must 
be promoted to satisfy non- vacuum equations of motion. 
The resulting theory will be a bi-scalar tensor theory of 
gravity, with the two scalars given by tp^ and ■0^. 



B. Bulk Geometry 

It is instructive to know the behaviour of the scalar 
field ip and warp factor w in the bulk. For instance, a 
relevant question is what the conditions are for a singu- 
larity to be present in the bulk? Let us start by noting 



that equations H43|) and (|44(l give us the way in which the 
scalar field tp and the warp factor uj behave as functions 
of the proper distance, d(^ ~ Ndz, in the bulk space. 
If Ub > 0, then the warp factor will decrease in the z 
direction, while if Ub < then the warp factor will be 
increasing. Similarly, if O^Ub > 0, then the scalar field 
will increase in the z direction, and if B^Ub < 0, then 
it will decrease. Moreover, if Ub{iP^) > then the first 
brane has a positive tension while if UB{tp^) < then 
it has a negative tension. The situation for the second 
brane is similar: if C/_b('0^) > then it has a negative 
tension while if Usi'^p'^) < then it has a positive ten- 
sion. In this way, it is possible to see that in general the 
two branes will not necessarily have opposite tensions. 
Figure 121 sketches the different possible behaviours for ip 
and to as functions of z. 



(b) 



V 




(c) 

4^ 



(d) 




\ 




\ 


^^'^ 


/ m 


\ 




> 



FIG. 2: The figure shows the four different possibilities for 
the behaviour of tp and w as a function of z, according to the 
signs of d^ Ub and Ub ■ The parameters of the figures are as 
follows: (a) a^(7s > and Ub > {01 a > 0). (b) O^Ub > 
and (7s < (or Q < 0). (c) O^Ub < and Ub > {or a < 0). 
(d) di,UB < and t/s < (or a > 0). 



From equation H44() we see that the infinitesimal proper 
distance, d(^, in the extra-dimensional direction can be 
written as: 



dC 



OUb 

dip 



dtp. 



(49) 



Observe from this last relation that extremum points of 
the potential Ub, given by the condition D^Ub — 0, 
can be only reached at an infinite distance away in the 
bulk. This is easily seen from the fact that dC -^ +00 
as dtiiUB — *■ 0. (The asymptotic behaviours at the ex- 
tremum points are, w — > if Ub > and w — > -|-oo if 
Ub < 0). Thus O^Ub cannot change sign in the bulk 
space and, therefore, tp{z) is a nionotonic function of z. 
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This is not the case, however, for the warp factor w. This 
can increase or decrease according to the sign of Ub(iP) 
for a given value of i/j. The fact that ipi^) is monotonic 
in the bulk space allows the possibility of parameterizing 
the bulk with ^ instead of z. This is an important result 
that will be heavily exploited in the rest of this paper. 

Finally, let us examine the possibility of having singu- 
larities at a finite position in the bulk. From equation 
(|47|l we can see that singularities will appear in the bulk 
whenever uj{z) = or uj{z) = +00. In general, these 
singularities will be located at points, Zqo, where the fol- 
lowing integral diverges: 






1/.1 



(50) 



Let us designate ip^c = '>p{zoc) to the value of the scalar 
field at which H5U|I becomes divergent. Then, from equa- 
tion H49|l. it is possible to see that for a singularity to 
be located at a finite proper-distance in the bulk, the 
following additional condition needs to be satisfied: 



Vo 



^1 



OUb 

dip 



dip 



< -1-00. 



(51) 



Now, notice that the integral in (|50|l will diverge either 
if d^UB{tpoo) -^ or UB{ipoa) -^ ±00. We have already 
studied the first case, which necessarily happens at in- 
finity. The second case, Ub{iP 00) ^ ±00, will be asso- 
ciated with singularities of the type uj —^ 0. For exam- 
ple, the only way to approach a singularity of the type 
Ub{'^) —^ ±00, with z — > z^, that is from the left, will 
be with O^Ub > 0. This corresponds to UBi^P) -^ -foo, 
and therefore to a warp factor going as w — > (the same 



argument follows for z 



3). This means that the only 



type of singularity, at a finite location in the bulk, will 
be of the type cu ^ 0. 

Summarizing, we have seen how the bulk geometry of 
the vacuum static system is completely determined by 
the behaviours of ip and u. -0 is a monotonic function of 
z, while w can increase or decrease depending on the sign 
of Ub- The only singularities possible in the bulk are of 
the form a; ^ 0. In general, this is a good reason to con- 
sider two-brane models in the context of BPS-systems: 
the singularity can be shielded from the first brane, and 
made to disappear from the bulk, with the presence of 
the second brane. However, the second brane can be at- 
tracted towards the singularity and eventually hit it. 



C. A Few Examples 



constant. As already said in Section Q] this form of the 
potential is motivated by heterotic M-theory, and is also 
predicted in five-dimensional supergravity. Let us first 
consider the gauge A^ > a constant, and, for simplicity, 
choose the positions of the branes at zi = and Z2 = r, 
with r > 0. Then, it follows that the solution to the 
system is given by: 



^P{z) = \n\~{Nz + s)a^Vo], (52) 

a 

w(z) = [-(iVz + s)aVo]*^exp(V'V4a), (53) 



where N and s can be expressed as: 



N = -- 






-a'lp' 



^-Qi/i^ 



(54) 
(55) 



For the singularities, there are two generic cases. When 
Vo > 0, the first brane is a positive tension brane and the 
second one a negative tension brane, and a singularity of 
the type w — +00 will exist at z = —00, while a singu- 
larity of the type w = will exist at the finite position 
Zoo = -s/N = r e-°"l'' / {e-°'l'' - e""'''") > r. Thus, it 
is not difficult to see that if atp'^ — > +00, then the sec- 
ond brane hits the singularity uj — 0. Figure |21 shows the 
solutions for 4'{z) and uj{z) in the present case (with the 
choice a < 0). In the case Vb < 0, the first brane corre- 
sponds to a negative tension brane and the second brane 
to a positive tension brane, with singularities happening 
at the same coordinates but with opposite signs. 



Singularity 




FIG. 3: The figure shows the solutions of tp{z) and uj{z) for 
the exponential case, when the gauge A'' constant is adopted. 



To gain some experience with BPS-configurations, let 
us briefly study the static vacuum solution for a few 
choices of the potential Ub- In particular, it will be useful 
to see how the Randall-Sundrum model arises as a par- 
ticular case of the present type of system. Let us start by 
analyzing the case Ub = Vqc"'^ , where a is an arbitrary 



The Randall-Sundrum scenario is easily obtained by 
letting a — !■ 0. In this limit, the above solutions can be 
reexpressed as: 



u;{z) 



0, 



-fiNz 



(56) 
(57) 



where /x = {1/4:)Ub- In this case, the singularities are 
located at infinity. Note that in this case the only degree 
of freedom is A^, the radion field. 



IV. LOW ENERGY REGIME EXPANSION 

We are now in a position to deduce the equations of 
motion governing the low energy regime. They consist in 
a linear expansion of the fields about the static vacuum 
solution found in the last section. As we shall see, these 
equations can be put in an integral form, allowing the 
construction of a systematic scheme to obtain solutions, 
order by order. The zeroth order solution of this ex- 
pansion corresponds to the static vacuum solutions with 
the integration constants ip^ and ip'^, and the metric g^^, 
promoted to be space-time dependent. 



A. Lo'w Energy Regime Equations 



the potentials. Additionally, 5^^ depends only on the 
space-time coordinate x. Note that if there is no matter 
content in the branes and u — vi — V2 =0, then we 
can take /i^jy — ip — (f> — 0, and the fields g^^, "0 and N 
of the previous definition would correspond to the static 
vacuum solution discussed in the last section, with iJjq 
and tJjq the two constant degrees of freedom. Therefore, 
the functions h^^,, ip and (f) are linear deviations from the 
vacuum solution of the system. Since we are interested in 
the low energy regime, we consider the case |/i^y| <^ |ff/ji/|, 
\f\ ^ iV'ol and \<j)\ <C 1. Now, if we insert these previ- 
ous definitions in the equations of motion 113|l . H15|l and 
(|19|l . and neglect second order quantities in terms of h^^, 
ip and (j) 'W6 obtain the equations of motion for the low 
energy regime. First, equation (|13|l leads to: 

dUf 



iVot/B(V'o)rX/3 + 2iVo— ^^' 



-m 



2Un 



dipo 



9^0 



'P 



-^X. (65) 



Let us develop the equations for the low energy regime. 
To start, assume that V'o, ^o and A^o are functions of x 
and z that satisfy the BPS-conditions: 



1 






No 



NoUsi^o) 
OUb 



dipo ' 



(58) 
(59) 



and let us define the bulk scalar field boundary values 
as V'o — "00(^1) and -00 = ■^0(2:2) (note that they are 
functions of the space-time coordinate x) . The solution 
for the warp factor, luq, is found to be: 



Equation 
tions: 



51 leads to the 4-dimensional Einstein's equa- 



K. 



NoUb{^o) 



dNo dUE 



di)Q d^po 



-7VoC/b(V'o)0'5m- 



3 -^ dUB , . 



--AT^ 



2f/n 



dtpo 



f 



/V2 

g,.. = 2^y^,. (66) 



OJn 



And finally, equation H19|) [with the help of equation l|13|) ] 
leads to: 



LOo{z,x) 



exp 



■4>o 



a 



K 



-\^)d^ 



(60) 



where the z and x dependence enter through the func- 
tions ipQ^x) and tpoizyx) in the integration limits. It will 
be convenient to define the warp factor between the two 
branes as w = iOQ{z2), or more explicitly, as: 



uj{x) = exp 



1 



a 



< 



-1(0) dV' 



(61) 



Now, we would like to study the perturbed system about 
the static vacuum solution. With this purpose, let us 
define the following set of variables, h^^, V and 0, as: 




N 






+ hf^i,) 



liv I , 



(62) 
(63) 
(64) 



where g^^, ip and N satisfy the equations of motion (|13|l . 
(|15|l and (|19|) . taking into account the presence of mat- 
ter in the branes and supersymmetry breaking terms in 



^ 



-Nn 



NoUb{^o) 



dNo OUe 



dtpo dtpo 



^ 



d^po 



N'o 



dijjo d% 



NS 



2 aV'o ^ 

1 



aoX . (67) 



In the previous equations, Uq = U{^o). Equations H65(l 
and (|66|l correspond to the linearized Einstein's equa- 
tions, while equation H67|l corresponds to the linearized 
bulk scalar field equation. In the previous set of equa- 
tions we have defined, for the sake of notation, the func- 
tions X, y^^ and 2, in the following way: 



X = uo' 



{di,f 



4 

-R + u 



y^v = G^,^,^—{gf,,UN-^^'Q,N) 



Z = UJn 



'^^-^^^^^^(^^^'^) 



(68) 

(69) 
(70) 
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X = 


= Xo 


3^p. = 




z -- 


= ^0- 



Note that in the functions X, y^j^^, and Z, the quanti- 
ties g^iy, "0 and iV are not exphcitly expanded. Their 
expansions can be considered in the following way: 

Z, (71) 

where Xq, Y^^^, and Zq are the zeroth order terms in the 
expansions of A", 3^^i. and Z (that is, they are constructed 
from ^^^, il)Q and Nq) while X, Y^^i, and Z are the terms 
which contain linear contributions from fe^i^, tp and 4>. 
Their specific form are given in Appendix^ We can now 
replace the new set of variables in the Israel matching 
conditions and the scalar field boundary condition. At 
the first brane, cti, these take the form: 



h'. 



3 



gt^i^T'^h'^f} = i^NoUB{ipa)<j)9tiu 



2 d^o 



'^Ub ~ ,3 . 2 AT T. 1 



(72) 



and 



^' = iVof;^0 + iVo^^ + iVo|f, (73) 

dipo d% dipQ 

where, quantities like ipQ and A^O: not explicitly evaluated 
at the boundary, must be evaluated a,t z = zi. Mean- 
while, at the second brane, 172, the matching conditions 
are: 

3 
/^U ~ 9^i,^9°''^K0 = 7;XoUb{i1^o) (j^gf,^ 



2 (9^0 



V3a"> + 2^QV2gt,i' + K^NoT^^uj„ , (74) 



and 









dv2 



(75) 



where again, quantities like ipo and A^o rnust be evaluated 
at z = Z2 . Note that in the boundary conditions the new 
set of variables are linearly proportional to the energy- 
momentum tensor at the branes. In other words, de- 
viations of the boundary conditions from the BPS-state 
generate the existence of fields h/^n, , (p and (j> in the bulk, 
as expected in the low energy regime. In order to solve 
the equations of motion it is useful to rewrite them in the 
form of a set of first order differential equations in terms 
of the dz derivative. This can be done by defining new 
fields ^^i, and x to be sources of the fields h^^, (p and </>, 
in the following equations: 



K. - ~g,^r^h'^p - :,NoUb{^o) c^-g,.. 



Nn 



OUb 
dipo 



3 -^ dUB . 
'2'^°9V^^^- 

' - A'O ^ ,y (f 



dlP'^ 



e^-, (76) 
X- (77) 



With these definitions oi ^^i, and X: it is found that equa- 
tions H()5|l . Ht)6|) and (|67|l can be reexpressed in a much 
simpler way. Respectively, these equations are found to 
be: 



d-00 



""^^ 



1 



A^ot/s(V'o)e 



— 2"^' 

^0 



'u;f 



No 






2NoLJ^y^., 



1 



aoAot^o ( Z - -aoX 



(78) 
(79) 
(80) 



where ^ = g'^^S.^v ■ Note that the left hand sides of equa- 
tions (|79|l and H80|l consist of first order derivatives of ^^jy 
and X in terms of z, while the right hand sides depends 
on second order derivatives in terms of the space-time 
coordinate a;. The boundary values of £,^1, and x, at the 
positions of the branes, can be computed with the help 
of the boundary conditions H72(l - H75|l . These are: 



eM^(^l) 


= -^Novig^,^ ~ kINqT^^, 


(81) 


X{zi) 


9-00 


(82) 


at position zi, and 






^^^^(^2) = 


-NoV2g^, + 4NoT;^,CJ-^, 


(83) 


X(^2) = 




(84) 



at position Z2. (Recall that uJo{zi) — 1 and 0^0(^2) — 
uj). We should not forget at this point the additional 
equations for the matter content at both branes. They 
come from equation (|4U|I and, with the present notation, 
are given by: 



VT^ = 0, 






{^,u^)r^Ti^ 



(85) 
2{V''Q)TlJ\, (86) 



where covariant derivatives V^ are constructed from the 
(jfj^i, metric. Additionally, it is worth mentioning that the 
projected Weyl tensor E'^^ and the loss parameter Aip 
can be written in terms of the linear variables as: 



E,, 






2No 

3, OUb X 






^fJ.l' offM" 



~a/3 



^^M^"" - 






lg^,{d^f + ±W^W,N 



(87) 



A0 = — 

No 



No 
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One of the main features of equations H79(l and (|80|l is 
that now they can be put in an integral form. That is, 
we can write: 



u, 



3. 



j^Ct^Az) = -^Qt^i^vi - kIT^^ 



+2 / dzNoUJ^y,,,, 



(89) 






+ / dzaoNouj^ ( Z - 2«oA' 1 , (90) 

where we have used the boundary conditions at the po- 
sition of the first brane. To further proceed we must im- 
plement a systematic expansion order by order in which 
the vacuum solution serves as the zeroth order solution. 
We develop this in the next subsection. 



B. Low Energy Regime Expansion 

Equations H89|) and (|90|) are two integral equations of 
the system, which suggest the possibility of introducing 
a systematic expansion about the vacuum solution, order 
by order. In the last subsection we saw that the matter 
content of the branes are sources for the linear deviations 
hf^^, (fi and (j) on the bulk. Therefore, it is sensible to 
study how the bulk is affected by the matter distribution 
on the brane, at different scales. In this way, let us con- 
sider the following expansion for the energy-momentum 
tensor of the matter content on the branes, as well as for 
the supersymmetry breaking potentials: 



and 



7.1 


- T^ (0) , 7.1 (1) I 7.1 (2) _| 


(91) 


Vl 


= .(°)-K.w+.f) + ..., 


(92) 


7.2 


_ 7.2(0) ,7.2(1) ,7.2(2) _l 


(93) 


V2 


= .f +««+.f + .... 


(94) 



The parameter of the expansion is dictated by the scale at 
which each term becomes relevant in the physical prob- 
lem of interest. Naturally, the expansion above will in- 
duce an expansion of the source functions ^p,y and Xj 
given by: 



f _ di) , f(2) , ... 
X = X(i) +X(2) H , 



(95) 
(96) 



in this way, the boundary conditions (|81|I - H84|) . can now 
be written, order by order, as: 

^l:^'Hzi) - lNov<C^~g,^~KlNoT;^!^^\ (97) 



X{i+i)izi) = No 



(0 






(98) 



at position zi, and 



ij:^'Hz2) = -Novi'^~g,, + 4NoT^J^co 



X(i+i)(^2) = N^ 



dv. 



(i) 



dipo ' 



(99) 



(100) 



at position Z2. Note the convention whereby the indicies 
of the left hand side are raised by one unit in relation 
those on the right hand side. Continuing with the con- 
struction, we must also consider the expansion of /i^^, tp 
and 6: 



''^U 

V 






"'fJ.U 



(101) 

"^(i) + '^(2) + • ■ ■ , (102) 

(j) = (/)(i) -K (/)(2) + • • • . (103) 

They are defined to satisfy the following first order dif- 
ferential equations with sources £,i^l and X(i) ■ 






<pU 






3 - , dUB 



Nn 



dUu 

dipo 



(0 



(9Vo 

d^UB 



^(») 



flu ' 



N, 



d^'o 



— V(t) =X(i)- 



(104) 
(105) 



Finally, recall that the quantities X , y^n and Z defined in 
equations H71|l depend on the fields /i^^, (p and (f>. There- 
fore, we must consider the following expansion 









X(2)- 



Zq + Z(i) + Z{2) + 



(106) 

(107) 
(108) 



where the index "0" denotes the dependence on the ze- 
roth order quantities, A^Oi V'o and wq, and the index i > 1 
denotes quantities that depend linearly on /i}j,y, Lpu\ and 
(\)ii\ . The precise form of the expanded functions X( j) , 

y^J; and Zu\ are shown in Appendix IXI Observe that it 
follows from equations H85|l and H86|l that the i-th order 
term, T^i/ , in the expansion of the energy momentum 
tensor for the matter living in the brane, satisfies the 
following conservation relation: 



^1/7.1(1) Q 



Y7 1/7.2(1) 



^\iy^r'T^ 



r.2(j) 



(109) 



(110) 



-2(V''c^)T^^,(^ 

With all these previous definitions we can now cast the 
equations of motion in the following way: 

X(j+i) - q'^s(V'o)4(j+i) = — 



(9Vo 



(111) 



f^;^(»+l) 



2iVoc.2r«, (112) 



«o^X(»+i) 



QoA^Ot^O ( Z{i) - -QoX(i) 



(113) 
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This last set of equations shows the desired low energy 
regime expansion. It states that ^)jV and X{i+i) can be 
solved in terms of the lower order quantities X/i\ , yj^^ and 
Z(j). At the same time, the functions A'(i), 3^^*j and Z(j) 

can be solved in terms of S^^Jll and xu) j '<^^ evident from 
their definitions. Note that the former means that we can 
compute E^^ and A^ to any desired order, as indicated 
by equations lP?7JI and (jHHJ- In this way, starting with 
the zeroth order solutions ipo, Nq and g^^, we can arrive 
at any desired order for the full solutions to the bulk 
equations ip, N and g^^. The present method is similar to 
the one introduced by Kanno and Soda for the Randall- 
Sundrum model [23) 1^ as well as to other schemes |^ 
139 . liol liH for dilatonic brane- worlds. 

We mentioned that the expansion parameter was dic- 
tated by the scale at which each term of the energy- 
momentum tensor expansion was relevant for the physical 
process of interest. In terms of the bulk quantities, one 
finds that the effect on the variation of the scale in the 
extra-dimensional direction is related to the variation of 
the scales in the space-time direction as follows: 



h" 



7V2 

-§-n/i(,), 

^0 



(114) 



That is, second derivatives of the (z+l)-th order variables 
in terms of z are of the same order as the second order 
derivatives of the z-th order variables in terms of space- 
time coordinates. 

As in the previous case, we can rewrite equations (|112|) 
and Ijll^l) in an integral form: 






No^^' 



U), 







+ 



2 dzNoLulY^, 

dvi 
Wo 



ao^rpX(»+i) = ai 



/ dzao 

J Zl 



NqUJq ( Z(i) - -Q;oA(i) 



(115) 



(116) 



This form of the expanded equations of motion are highly 
relevant; they are the basis for the derivation of the ef- 
fective theory developed in the next section. 

Summarizing, in this section we have developed an ex- 
pansion procedure to solve the entire theory in a system- 
atic form. As we said, to solve this system of equations 
we need to know the form of the zeroth order moduli 
fields ipQ and f/^^, and the metric g^j/ and, therefore, an 
effective theory for these fields is required. In the fol- 
lowing section, we are going to deduce the 4-dimensional 
equations of motion governing these functions, as well as 
those for higher order terms in the expansion. 



V. FOUR DIMENSIONAL EFFECTIVE 
THEORY 

In the previous section we have shown how to construct 
the low energy regime as a consistent expansion. In this 
section we see how to define the four-dimensional effective 
theory at the position of the branes starting from the 
expansion above. Concordant with the expansion, the 
effective theory must be defined order by order. 



A. General Case 

The effective equations governing the variables at the i- 
th order can be obtained by evaluating the integral equa- 
tions H115|l and H116|) at z = Z2. In other words, the z-th 
order equations are given by: 



dzA^o-o%':) = ^[^V^-<^]5.. 



-^? 



7^2(0(7,2 ,7^1 (») 



(117) 



dz aoA^owg ( Z(i) - -a^X^i-^ \ = ol2L 



I dv. 



{i) 



di>l 



-ai-fr, (118) 

Additionally, another equation can be obtained by eval- 
uating equation (|lll|l at the boundary position zi. The 
resulting equation is: 



A(,)(zi) = 2 



(') 



dUB dvl 

di>l d-ipl 



2UB{i^l)vf^ 



^^usi^Drx^'^- 



(119) 



Note that a similar equation would have been obtained 
by evaluating (|lll|l at the boundary position Z2, however, 
this would not give a new equation of motion. Equations 
(|117|l - (|119|l are the desired equations of the 4-dimensional 
effective theory. Solving these equations for the z-th order 
variables allows us to obtain the effective equations for 
the {i + l)-th order variables, after correctly integrating 
(flMjl and ifTUsll . The form in which equations ((TTTjI - lfTT^ 
are presented is, at this level, abstract and it is difficult 
to appreciate the effective theory in more familiar terms. 
In the next subsection we analyze in detail the effective 
theory at the zeroth order in the expansion (i = 0). 



B. Zeroth Order Effective Theory 

It will be particularly interesting to analyze the form of 
the effective equations at the zeroth order in more detail. 
In this case the metric at the first and second branes are 
conformally related. In terms of the treatment developed 
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in the previous section, the effective theory is: 



,2^0 3_ 

^2 r 



dz A^o^o^u^ = -;9piu 



-4 (0) 

Ld V2 — V\ 



T2(0)-2 , ^l(O) 



I dz aoA^o^o ( -^0 - 2*^0^0 ] = ^2^ 






dv{ 



(0) 



^o(^i) 



dUsOvf^ 1 (0) 



(120) 



(121) 



' di'l av-o 



-C/s(V'^)5''''Tiio). (122) 



In the foUowing, we shall omit the "0" index denoting 
the zeroth order terms. To rewrite this theory in more 
familiar terms, we need to integrate the terms on the left 
hand side of equations (|12(J|I and p21|) . In order to do 
this we need to evaluate the terms at the boundaries. 
Hence it is particularly useful to express the theory in 
terms of the moduli fields -0^ and ijP' . It is then possible to 
obtain the following effective theory f Appendix 151 shows 
how to derive the next results). The Einstein's equations, 
obtained from (|120|) . are found to be: 



.M [^2j,, 



2 






-^9,^V, (123) 



while the moduli fields equations, obtained from equa- 
tions I12U|) and (|122|) . are found to be: 



D^" = --g^'^^id^^^d.r + 11"'^, 



;7 












(124) 



In the previous expressions, the index a labels the posi- 
tions 1 and 2. Additionally, note the presence of a con- 
formal factor n^ in front of the Einstein's tensor G^^, 
given by: 



n' 



ftp' 
k / diij 



OUb 

dtp 



(125) 



where uj is given by equation (|47|l . The coefficient k is 
an arbitrary positive constant with dimensions of inverse 
length, which has been incorporated to make fl^ dimen- 
sionless. The symmetric matrix jab is a function of the 
moduli fields, that can be regarded as the metric of the 
space spanned by the moduli in a sigma model approach, 
with 7'^'' its inverse. The elements of jab are given by: 



711 = "1 ■ 



Uei^Jj^] 



2 



(126) 



722 



7l2 



iu^k 



Ub{4'^) 



2V 



—a, tto 



C/b(02)' 



(127) 
(128) 



with 721 = 712. Associated with the above metric, we 
have defined a set of connections 7^^. These are given 
by: 



1 



Ibc 



t7 



ad 



djbd , djdc djbc 



+ 



dip" dip'' dip" 



(129) 



Finally, we have also defined an effective potential V 
which depends linearly on the supersymmetry breaking 
potentials u, vi and V2. This is defined as: 



V = 



dip 



fdUE 



l/;l 



\ dip 



4 

U! U - 



2fc 



LO V2 — Vl\ 



(130) 



Finally, we can not forget the matter conservation rela- 
tions. These take the form: 



VT^ = 0, 



(131) 



V%'. = ^[{V,u;)r^Tlp-'2{V''u^)Tl,]. (132) 

The fact that the energy-momentum tensor of the second 
brane does not respect the standard matter conservation 
relation is due to the frame chosen to describe the system. 
The generic form of the theory displayed by equations 
(|123|l and (|124|l is of a bi-scalar tensor theory of gravity, 
with the two scalar degrees given by tp"^ and tp"^ . The 
above set of equations can be obtained from the following 
action: 



S = J^J d'-^V ~9 



^''R~-^r''jabd^,rd.tp''-\v 



-5i[*i,.9^^]+52[*2,c:'^5/.^], 



(133) 



Equation H133() is an important result. It can be shown 
that the effective action H133|l corresponds to the moduli- 
space approximation 42] , where the relevant fields of the 
theory are just simply the free degrees of freedom of the 
vacuum theory, promoted to be space-time dependent 
fields. A relevant aspect of this theory, is that the equa- 
tion of motion for a moduli field -0° will depend linearly 
on the trace of the energy-momentum tensor of the brane 
a but not on the one belonging to the opposite brane. 
More precisely, it can be shown that p24|) has the form: 



Uip'^ = 



kK5 dUs y^ 1 
6 9V^ ^"^ 

fcKs dUs 
6 dtp'- 



-u g 



2~fll'rp2 



(134) 
(135) 



This means that, in the low energy regime, the moduli 
fields are driven by the matter content of the branes (re- 
call that the moduli are parameterizing the positions of 
the branes). This behaviour of the moduli fields is inde- 
pendent of the frame choice (see next subsection) . 
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C. Einstein's Frame 

Note that in equation H133() the Newton's constant de- 
pends on the moduli fields. This theory can be rewritten 
in the Einstein frame where the Newton's constant is 
independent of the moduli. Considering the conformal 
transformation: 



9fj.v *■ 9 111' — ^' 9iJ.vj 
we are then left with the following action 



S = 



1 



kh 



dS 



-9 



R 



3 



9''''iabd^rdur 



(136) 



:V 



+Si[^^,Alg^,]+S2[^2,Alg^,], (137) 

where now the sigma model metric 7a& is given by: 



711 



2a: 



722 = 2q, 



t/l(^l) 

-2 '^ ^2 



7l2 



2aia2 



UB{i'')A^' 



l + :^C^i3(^')V 



2fc 

1 

2fc 



i-242 /|2 



C/s(^i)t/s(^2)- 



(138) 
(139) 
(140) 



As usual, 7"'' is the inverse of 7af,, and j^b ^^^ ^^^ connec- 
tions. It is possible to show that ^ab is a positive definite 
metric. To work in this frame is particularly useful and 
simple. Additionally, we have defined the quantities Ai 
and A2 (which are functions of the moduli) to be: 



AI' 


= n' 










J 4,^ 


V 5^ . 


-1 
jexp 


1 
~2 


A- 


= n^Lu-^ 










= k d^ 


fdUB-^ 


-1 
jexp 


1 
~2 



a~ di/j 



a~^d'tp 



(141) 



(142) 



Also, the potential V is now found to be: 

fV" /orr S -1 



V 



-n- 



dip 



dip 
-2k [Ajv2 - Ajvi] 



1/ji 



4 

U! U 



(143) 



The present form of the theory, can be further worked 
out to put the sigma model in a diagonal form with the 
redefinition of the moduli fields. However, in the present 
work we shall continue with the theory in the current 
notation. Let us finish this section by expressing the 
equations of motion in the Einstein frame. The Einstein's 
equations are: 



G^. = -^iab{]^9i..9''''d^rdp^l^'' - d^rdu^" 



(144) 



and the moduli fields' equations are: 



□V^' - -9'"'ll{d,i''){d.r) + U'^^^ 



2 ' dip" 



knl OUb 



6 dij: 



1^1 9 ^121'^ 



□^' = -9''-llid,^'')idur) + \l''§^ 



kni OUb 



Al g'^'T,'^ 



(145) 



(146) 



6 9V^ ^ M^- 

Additionally, the matter conservation relations now read: 

V%^. = -^[(V,Ai)5"^ri^-2(V''^i)r^i,], (147) 
V%'. = -^[(V,A2)5"^r2^-2(V''^2)r2J. (148) 

In the next sections, we are going to study this the- 
ory in more detail, analyzing a few examples (including 
the Randall-Sundrum case) and studying the cosmolog- 
ical evolution of the moduli for late cosmological eras. 
We shall pay special attention to the observational con- 
straints that can be imposed on this model. 



D. A Few Examples 

In this subsection we shall briefly analyze two well 
known cases, namely, the exponential case Ub{iP) = 
Vbe"'^, and the Randall-Sundrum case, which can be 
derived as particular cases of the formalism described 
above. 



1. Exponential Case 

This case is particularly simple, and has been studied 
in detail in several previous works [2J, [23, 123, l^l, |23, 

US 113. The potential to consider is t/B(^) = Vbe"'^, 
with a being a constant parameter of the theory. Here, 
it is possible (and convenient) to parameterize the theory 
using a new set of fields. Let us assume that Vb > and 
consider the following transformations: 

2 



„2*i 



cosh $9 



6^*1 sinh^ $2 



l + 2a2 

2 



^(l+2Q^)i/)V2a 



-(l+2ct)^^ /2a 



l + 2a2^ • (149) 

If Vb < then we should interchange ip'^ and ip"^ in the 
previous definition of <&! and $2- Inserting this notation 
back in the theory, it is possible to obtain the following 
effective action: 

Q ^ [ A^ r-\- 12a2 



knl 



R 



6 



;{d^2f 



l + 2a2 



:V 



{d^iY 



l-h2a2'" '' 4 
+Si[^i,Alg^,,]+S2[^2:Alg^,\ 



(150) 
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where now the coefficients Ai and A2 are functions of the 
fields $1 and $2: given by: 



isotropic Universe. We can derive the equations using 
the flat Friedmann-Robertson- Walker (FRW) metric: 



4? = 


_ g-2*,]^ 


'1 
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[2 


42 - 
^2 - 


- p-'^'iiYi 


'1 
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[^ 



(1 + 2a')e^*i cosh" $2 



(1 + 2a2)g2<i.i gjj^j^2 ^^ 



• (151) 



This form of the theory was already obtained in ;3Q | , in 
the moduli-space approximation approach, and extensive 
analysis to it have been made. 



2. Randall Sundrum Case 

Finally, it is also instructive to examine the equations 
for the Randall-Sundrum case. This can be obtained as a 
limiting case of the previous example, by letting a — > 0. 
The only field surviving the limiting process is $2, and 
the effective action is: 



1 



S ^ -^ I <^ ^\/^ 



R-6[d^2) 



^y 



+Sx{^^,A\g,,,\^S2{^2.J^29^.v\ 



(152) 



where now the coefficients A\ and A2 are functions of the 
field $2 alone, and are given by: 



A\ = ^cosh^$„ 



Al 



2k 
2k 



sinh $2. 



(153) 



In this case, however, the field is not directly coupled to 
matter. 



ds^ = g^^dx^'dx" = -dr + a^t)d^jdx'dx^ . (154) 

Here, g^^i, corresponds to the Einstein's frame metric, 
conformally related to the physical metric, g^^i,, of the 
brane frame. Using the metric H154() . jointly with the ef- 
fective action H137(l . we are left with the following Fried- 
mann equations: 

" -'^At{p,+3p,)~^-^Ai{p2+3p2) 



a 12 



12 



-^^abr^' + lv, 



6 6 8 



rK 



(155) 



(156) 



where V is given in equation (|143|) and H = a/a is the 
Hubble parameter. The moduli fields equations, on the 
other hand, are given by: 



IJ.aJ.b KlcdV 



^i+3H^i = -^l^^-^b__^ ^^^ 



-ifc«^^At(pi-3pi), (157) 



^' + 3H^' = -^a^ri^'-l^'^^ 



ifc.29C/B ,4 



knl^A^2{p2-ip2). (158) 



Additionally, the matter conservation relations can be 
read as: 



VI. LATE-TIME COSMOLOGY 

A simple and important application of the formalism 
developed above is the study of cosmological solutions 
for the most recent epoch of the Universe. In simple 
terms, current observations reveal that the moduli fields 
are not relevant for the present evolution of the universe 
and their dynamics is strongly suppressed. In this section 
we shall analyze the possibility of stabilizing the moduli 
fields V'^ and ijP' . Interestingly, it will be found that, in 
order to stabilize the system, constraints must be placed 
on the global configurations of branes, particularly, on 
the class of possible potentials and also on the position 
of the branes with respect to the background. 



A. Cosmological Equations 

Let us derive the cosmological equations of motion of 
the present system, valid for a flat, homogeneous and 



pi+3i/(pi+pi) = -3-i(pi+pi), (159) 

P2+iH{p2+P2) = -3-^(p2+P2). (160) 

A2 



The above set of equations are the corresponding equa- 
tions of motion describing the cosmological behaviour of 
the brane system. Since the sigma-model metric is pos- 
itive definite, the moduli fields are not going to result 
in an accelerating universe. Thus, in the present con- 
text, the only way to obtain an accelerating universe is 
to consider supersymmetry breaking potentials. A more 
detailed analysis of this system, taking into account the 
full dependence on the extra kinetic terms coming from 
the sigma model formulation is complicated, and shall be 
omitted in the present work. However, since we are inter- 
ested mostly in the low energy regime of branes, we shall 
study the case in which the moduli are slowly evolving 
compared to the Hubble parameter. We examine this in 
the next subsection. 
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B. Late-Time Cosmology 



Observational Constraints 



Let us consider the case in which the moduH fields 
are slowly evolving compared with the Hubble parame- 
ter, that is \ip^\, IVj^I <C H, and that the supersymmetry 
breaking potentials are such that m, ui, W2 = 0. Then, the 
Friedmann equations of the system are: 



- = -'^At{p,+3p^)-^Atip, + 3p,), (161) 



12 



T2" 



2 kn^ 4 KKs 4 

H = ——A^pi + ——A2P2- 



6 



6 



(162) 



The moduli fields equations, on the other hand, are given 
by: 



V^l+Si/T/'l 



6 



dtp 



^^t(pi-3pi), (163) 



And the matter conservation relations can now be written 
as: 



Pi+iH{pi+pi) = 0, 

P2+iH{p2+P2) = 0. 



(165) 
(166) 



It is possible to appreciate that the moduli fields are 
driven by the matter content of the branes, while the evo- 
lution of the Universe remains unaffected by the moduli. 
In particular, if the universe is matter dominated, then 
there is an attractor for the moduli towards the extremes 
of the supersymmetric potential Ub ■ To be more precise, 
ip^ will be attracted to the minimum of Ub, while -0^ will 
be attracted to the maximum of Ub- This is sketched 
in Figure^ It is therefore sensible to assume that the 




> ¥ 



FIG. 4: The figure shows the case when the branes are mat- 
ter dominated. The moduli fields ij)^ and 'tp'^ are driven to the 
minimum and maximum of the superpotential Ub, respec- 
tively. 



present state of the Universe is very close to this con- 
figuration, where the potential is being extremize by the 
moduli. In the next subsection we assume that this is 
the case in order to place observational constraints on 
the model. 



We now turn to the task of placing observational 
bounds on the theory. In the last subsection we observed 
that the branes are cosmologically driven by the mat- 
ter content in them. In particular, the evolution of the 
moduli is such that they extremize the value of the su- 
persymmetric potential Ub- Therefore we expect that at 
late cosmological eras, such as the present era, the branes 
are far apart and their tensions nearly at their extreme 
values. This corresponds to an interesting global condi- 
tion on the configuration of the system, however, it can 
be constrained by current observations. More precisely, 
the measurement of the Post Newtonian Eddington coef- 
ficient 7 is constrained by the Very Long Baseline Inter- 
ferometry measurements of the deflection of radio waves 
by the Sun to be 3.8 10"* > 1 - 7 > -2.6 lO"'^ at 68% 
confidence level [43Ll44| . The parameter 7 is defined as: 



1-7 



?7" 



I + 11' 



2rf 



where 77^ — '^"''' rjarjb, and r/a is defined as: 



Va = 



1 dAi 

A^dtp'"' 



(167) 



(168) 



Observe that r]a is constructed from Ai, which is the 
factor that appears in the action for the matter content 
at the (Ti brane. This is because we are assuming, without 
loss of generality, that we are performing measurements 
from the cti brane. The quantity rj^ can be computed 
exactly: 



1 



7] 



I-^Ub^A-' 



Then, we finally obtain: 

2 ' 
3 



1 



7 



1 



2fc 



Ub{^j')A^' 



(169) 



(170) 



(Note that this result is independent of k). Recall that 
Ai is related to an integral over the entire bulk. This 
is a very important result: observational measurements 
constrain the global configuration of the brane system. 
For example, we can compute the 7 parameter for the 
exponential potential. In this case, we obtain: 



7 



31- 
2 



2a2 
1 



3 1 -h 2a2 



,(1+2q2)(V'-V')/2o 



(171) 



which agrees with earlier computations [3(1] . If the branes 
are such that the moduli are extremizing the potential 
Ub, then we have a(V'^ — "0^) -^ -l-C!0, and the second 
term will disappear in (|171() . Then, we end up with a 
constraint on a, the parameter of the model: 



1 



■7 



3 IH- 2a2 ' 



(172) 
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which gives \a\ < 10~^. Interestingly, this bound does 
not affect the global configuration of the brane, but only 
the parameter of the theory. This is certainly a problem 
since theoretically viable values for a are usually of the 
order of unity. In the Randall-Sundrum case there is no 
potential, and the moduli can not be driven to any stable 
configuration. In this case, the constraint takes the form: 



1 



2 , 

7 ~ -tanh $2. 



(173) 



Our results are applicable to many other aspects of 
brane-worlds not considered in this paper. For exam- 
ple, the development of the theory at the first order in 
the perturbed metric would allow an ideal background to 
study the cosmological perturbations and the CMB pre- 
dictions [43, llal . Additionally, many results where the 
exponential case Ub oc e"'^ was considered can now be 
extended to the general case of an arbitrary potential. 



Thus we see that in the Randall-Sundrum case a stabi- 
lization mechanism for the radion field is necessary in 
order to agree with observations. 



VII. CONCLUSIONS 
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In this paper we have studied several aspects of the 
low energy regime of EPS brane-world models. We have 
developed a systematic procedure to obtain solutions to 
the full system of equations, consisting in a linear expan- 
sion of the fields about the static vacuum solution of the 
system. As a result, it was shown that five dimensional 
solutions can be obtained at any desired order in the ex- 
pansion. In particular, the projected Weyl tensor .E^t^ 
and the loss parameter Aip can be computed with any 
desired accuracy. Additionally, and probably more im- 
portant, an effective 4-dimensional system of equations 
have been obtained. Concordant with the method, these 
4-dimensional effective equations must be considered up 
to the desired order in the expansion. For instance, we 
have analyzed in detail the zeroth order effective the- 
ory, which agrees with the moduli-space approximation 
p3 |. At this order, the metrics of both branes are con- 
formally related, and the complete theory corresponds to 
a bi-scalar tensor theory of gravity [equation (|137|l ]. The 
approach followed in this paper to obtain the effective 
theory is similar to the one followed in a previous work 
[3ll | , where the cosmological setup of the case Ub oc e"^ 
was investigated. However, the treatment in this article 
was more general in two senses: we did not restrict the 
form of the BPS-potential, and we deduced the complete 
Einstein's equations governing the system. 

We have also seen that the moduli fields -as defined 
in the zeroth order part of the theory- can be stabilized 
with the help of effective couplings between the mod- 
uli and matter, that arise naturally. This result comes 
from the fact that the equations of motion for the pro- 
jected bulk scalar field on the brane depend only on the 
energy-momentum of the respective brane. This allows 
satisfactory late cosmological configurations and enables 
us to place observational constraints on the model. An 
important result was the computation of the Post New- 
tonian Eddington coefficient 7 in terms of the moduli 
[equation (|170|l ]. This result is relevant since 7 is cur- 
rently the most constrained parameter of General Rela- 
tivity. For example, in the exponential case it was found 
that jaj < 10-2. 



APPENDIX A: DEFINITION OF X, Y^^ AND Z 

Here we define the zeroth order and the linear expan- 
sions of X, y^i, and Z of equations l(7T]l . The zeroth 
order quantities are: 

XQ^[diJof-^R + ^ujQ^nuuo+ujlu, (Al) 
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L^ui/ 



:9fii^idi>a)'^ - df.ipod^'ipo 
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39,.wo9^(A'"owo) 
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(A2) 
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■.g'"'V^(iVo9.7/.o) 
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(A3) 



Meanwhile, the linear terms in the expansions of X, y^^ 
and Z are given by the following expressions: 



X = 2 d^Jod^ - - (V^ V^/i'"' - Uh) 



VA/iQ/3)VpLi;o, 
1 



(A4) 



V = o (v<tV^/i;^ + v<,Vp/i;: - v^v^/i - n/i^^ 



2" 



3 

+g^„n(l) - V^V,.0 + g^„uja^g°''^daUJodp(j) 
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- 2~g^,r^Vo.hpx + ~g^..^xh)r^^pNQ, (A5) 

- 5^'V^0V,Vo - □¥> - 2^o"i.g''^9ACc'oap(^. (A6) 

To compute the i-th order terms of the hnear expansions, 
with i > 0, it is enough to add an index "i" to every hnear 
variable in the expressions above. 



APPENDIX B: COMPUTATION OF THE 
ZEROTH ORDER THEORY 

Here we indicate how to compute the integral in the 
left hand side of equation p2U|) . The integral present in 
equation H121|l can be solved in a similar way. First, it is 
important to note the following identity: 

4 
d^{Nf)UjQ) = -NoaoUJodf^^o - -— (9^wo)'- (Bl) 

Additionally, recall that we can parameterize the z coor- 
dinate using the monotonic zeroth order solution for the 
scalar field "00 (-2) • That is, we can write: 



, 1 fdUBY', 

'^^^ Tr\ '^rr ^^o- 



(B2) 



Using the previous equations, with boundaries ijjqIzi) — 
ip^ and ^"0(^2) — V'^i then the mentioned integral can be 



solved as: 



k I dzN^LolYll, 



f^'G^. + 5M-n«' - v,,v,o2 



^, 3fc /-^ (dUB\-\ 



(B3) 



where fi^ and ^ab are defined as in Section To finish, 
we should mention that in obtaining this result it was 
useful to notice that VP' does not only depend on -0^ 
and ip"^ through the integration limits in (|125|1 . but also 
through w^ present in the integrand. Recall that lu is 
normalized to be 1 at the position of the first brane, and 
therefore depends on ij)^ . This in turn means that the 
space-time derivative of il^ will have the form: 



d^n' 



1^2 _^_ 

2 Ub{V) 



_i LJ K 2 



5m^' 



(B4) 
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